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ABSTRACT 

We study the dependence of the galaxy contents within halos on the halo formation time using two 
galaxy formation models, one being a semianalytic model utilizing the halo assembly history from a 
high resolution TV-body simulation and the other being a smoothed particle hydrodynamics simulation 
including radiative cooling, star formation, and energy feedback from galactic winds. We confirm the 
finding by Gao et al. that at fixed mass, the clustering of halos depends on the halo formation time, 
especially for low-mass halos. This age dependence of halo clustering makes it desirable to study the 
correlation between the occupation of galaxies within halos and the halo age. We find that, in halos of 
fixed mass, the number of satellite galaxies has a strong dependence on halo age, with fewer satellites 
in older halos. The youngest one-third of the halos can have an order of magnitude more satellites 
than the oldest one-third. For central galaxies, in halos that form earlier, they tend to have more 
stars and thus appear to be more luminous, and the dependence of their luminosity on halo age is not 
as strong as that of stellar mass. The results can be understood through the star formation history in 
halos and the merging of satellites onto central galaxies. The age dependence of the galaxy contents 
within halos would constitute an important ingredient in a more accurate halo-based model of galaxy 
clustering. 

Subject headings: galaxies: formation — galaxies: halos 
mology: theory — dark matter 



large-scale structure of universe — cos- 



1. INTRODUCTION 

In the cold dark matter hierarchical model of struc- 
ture formation, the formation and evolution of galaxies 
are coupled with those of dark matter halos. Studying 
the distribution of galaxies inside halos can aid us in the 
understanding of the galaxy formation process and in the 
interpretation of clustering properties of galaxies. In this 
Letter, we investigate the dependence of the galaxy con- 
tents on the halo formation time using numerical galaxy 
formation models. 

Our study is closely related to recently developed 
models of galaxy clustering, namely, the framework of 
the halo occupation distribution (HOD) and the ap- 
proach of the conditional luminosity function (CLF). 
Both models are based on statistical descriptions of the 
relation between galaxies and dark matter halos. The 
HOD characterizes the relation in terms of the probabil- 
ity distribution P(N\M) that a halo of mass M hosts 
N galaxies of a given type, and the spati al and ve- 
locity distributions of galaxies within halos lling et al.| 
19981 iMa fc Frvl 120001 peacock fc Smith 1 1200(1 iSeliakl 
200(1 iScoccimarro et al.l200lHBerlind fc WeinberJ2002t 
Coorav fc Shetbll2002f) . Instead of P{N\M) for galax- 
ies of a given type, the CLF method uses the lumi- 
nosity function of galaxies within halos as a function 
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of halo mass to establis h the relation betw een galaxies 
and dark matter halos ijYang et alJ 120031) . HOD and 
CLF models have been applied to interpret t he observed 
galaxy clustering in a number of surveys (e.g.J Jing et a~ 
| 1998ll2002HBullock et alJl2002HMoustakas fc Somervil '' 
pnna lYangetaU 12003 Iva,n den Boscli i et a.U 
i Magliocchetti fc Porciariil 12003 lYan et all 12003 : iZheni 
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20041 IZehavi et al] 12004 12003 IQuchi et all I200E 
iLee et aTTl2005FICooravl2005ll , and with the help of these 
models our understanding of galaxy clustering has been 
greatly improved and refined. 

One key assumption in these halo-based models is that 
the statistical distribution of galaxies inside halos de- 
pends only on halo mass. This assumption has its theo- 
retica l origin in the excursion set formalism ijBond et alJ 
11991(1 . which predicts that a halo's assembly history de- 
pends only on its mass. It is inherent in any semianalytic 
galaxy formation model that uses halo merger trees based 
on the excursion set formalism. Using the high-reso lution 
Mi llennium iV-body s imulation, iGao et al] (|2005T ) (also 
see lHarker et al]l200"?)j) show that the clustering of halos 
at a given mass has a dependence on the halo forma- 
tion time, and the dependence becomes much stronger 
at halo masses <§; A/*, the nonlinear characteristic mass. 
Galaxy properties are expected to be correlated with the 
halo formation time, and to this extent, the environmen- 
tal dependence of the halo formation time would lead 
to a dependence of the statistical distribution of galax- 
ies inside halos on the environment as well as on halo 
mass. Once high accuracy modeling eventually becomes 
demanded by galaxy clustering data and/or by cosmo- 
logical applications, the current HOD and CLF models 
may have to be modified to take into account the envi- 
ronmental dependence. Given the age dependence of the 
halo clustering, the study of galaxy contents as a func- 
tion of halo age at fixed halo mass would shed light on 
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the extension of the HOD/CLF model. Moreover, the 
dependence of galaxy properties on the halo formation 
time is an interesting problem on its own, since it can 
aid us in our understanding of the galaxy formation pro- 
cess. 

In this Letter, we examine the dependence of the halo 
occupation properties of galaxies on the halo forma- 
tion time, using a semianalytic galaxy formation model 
(SAM) and a smoothed particle hydrodynamics model 
(HYD). In § 2, we introduce the simulations and galaxy 
formation models. We show a confirmation of the age 
dependence of halo clustering in § 3. The dependence 
of galaxy properties on the halo formation time is the 
content of § 4, where we present our results in terms of 
the CLF as a function of halo mass and formation time. 
In § 5, we give a summary and discussion. 

2. THE SIMULATIONS AND GALAXY FORMATION 
MODELS 

Two main simulations are used in our study. The 
one that t he SAM is ba s ed on is a P 3 M A-body sim- 
ulation in Ding fc Suto I l)2002|l with 512 3 particles in 
a cubic box of 100 /i _1 Mpc. A spatially-flat ACDM 
model is adopted with density parameters fl rn = 0.3 
and f^A = 0.7. We adopt the co ld dark matter power 
spectrum of iBardeen et alJ (|1986f> with a shape param- 
eter T = Q, m h = 0.2 and the amplitude <jg = 0.9. The 
softening length is 10/i _1 kpc, and the particle mass is 
M p ~ 6.2 x 10 8 /i _1 Ajf Q . We denote this simulation as 
SAM- ACDM. The gal axies of the SAM s ample are gen- 
erated by the model of Kan g et all (|2005Tl that explicitly 
follows the assembly history of each halo in the simula- 
tion. The model successfully explains a wide variety of 
observational facts, including luminosity functions in dif- 
ferent bands and the bimodal color distribution of galax- 
ies. 

The other simulation, denoted as HYD-Gadget2, is 
performed as smoothed pa rticle hydrodynamics with the 
TRE E-PM code Gadget2 ijSpringel 112003: iSpringel eTaD 
1200 following the evolution of 512 3 dark matter parti- 
cles and 512 3 gas particles in a cubic box of 100 /i _1 Mpc. 
It includes the physical processes of radiative cooling, 
star formation, and energy fe edback from galactic winds 
( Spri ngel fc Hernauistl I2003|L The cosmological and 
initial density fluctuation parameters of the model arc 
(fi m , Q A , n b , a s , n, h)=(0.268, 0.732, 0.044, 0.85, 1, 0.71), 
and the transfer function is computed from CMBFAST 
ijSeliak fc Zaldarriaga I ll99filL The softening length is 
also 10 /i _1 kpc, and the particle mass (dark matter +gas) 
is M p ~ 5.54 x 10 8 h~ 1 M Q . Galaxies are identified 
as gravitationally bound groups of stars and cold gas 
particles that are associated with a common local 
maximum in the baryon density. We note that the effect 
of dust extinction, which is included in the SAM model, 
is not taken into account in this model. 

In addition to the two main simulations, we also use 
a low-resolution simulation of 512 3 particles in a box of 
300 /i _1 Mpc only for the purpose of probing the age de- 
pendence of halo clustering at high halo mass. This sim- 
ulation (denoted as 300-ACDM) has the same cosmolog- 
ical parameters as SAM-ACDM, and the particle mass is 
M p ~ 1.67 x IO^/t^Mq. 

In all the simulations, ha los are identified with the 
friends-of-friends algorithm ijDavis et a □ with a 
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Fig. 1. — Halo bias factor as a function of formation time. The 
two integers in brackets indicate the range of particle numbers of 
halos (symbols have the same meaning in the two top panels) . The 
particle mass is labeled as M p in each panel. The open symbols 
connected by lines are the results in different age bins, while the 
filled symbols show the mean formation redshift and bias factor for 
all halos in the given mass bin. Poisson error bars are plotted for 
the two high-mass bins in the bottom panel. 



linking length 0.2 times the mean particle separation. 

3. THE DEFINITION OF FORMATION TIME AND THE 
AGE DEPENDENCE OF HALO CLUSTERING 

We define the halo formation time as the time when 
its most massive progenitor has exactly half of the final 
mass, interpolated between two adjacent outputs if nec- 
essary. The progenitor of the final halo is defined as a 
halo in an earlier output that has more than half of its 
particles found in the final halo. If the most massive 
progenitor in an output has more than twice the mass of 
the one in the earlier adjacent output, we excluded the 
final halo in our study since its formation time is difficult 
to measure because its main progenitor was likely to be 
bridge-connected to a massive halo. In the end, about 
10% of the halos are excluded. 

We divide halos in each halo mass bin into several(3, 
5, or 10) age bins with an equal number of halos in each 
age bin. The square of the halo bias factor is calculated 
by taking the ratio of the two-point correlation functions 
of halos and matter, averaged over scales 5/i _1 Mpc < 
r < 20ft,- 1 Mpc. 

The results of the dependence of halo clustering on 
the formation tim e are shown in Figure Q We con- 
firm the finding bv lGao et alJ ((2005) that old halos clus- 
ter more strongly than young ones. In particular, at 
~ 10 11 h~ 1 M ( 7 ) , the large-scale autocorrelation amplitude 
for the oldest 10% of the halos is more than 5 times that 
for the youngest 10% of the halos. For massive halo s, the 
dependence becomes weak. iWechsler et alJ l)2005l) find 
that the age dependence of the halo clustering reverses 
its trend for halos more massive than the nonlinear mass 
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Fig. 2. — Conditional luminosity function (CLF) and conditional stellar mass function (CSMF) as a function of halo age. Left two 
columns: CLFs from the SAM (left) and HYD (right) galaxy formation models. Results are shown in three mass bins indicated by the 
particle numbers [the particle mass in the SAM (HYD) model is M p =6.2(5.5)xl0 8 h _1 M Q ]. Central two columns: Same as the left 
two columns, but the CLF is separated into contributions from central galaxies (filled symbols) and satellites (open symbols). Right two 
columns: Same as the central two columns, but for the CSMF. The CSMF is not complete at low stellar mass, which causes the steep drop 
at the low stellar mass end. 



M* (10 13 and 4.6 x 10 12 /i^M© for the SAM-ACDM and 
HYD-Gadget2 simulations in our study), if halo concen- 
tration is used as an approximation of age. We do not 
detect such a clear trend, which could be buried in the 
large (Poisson) uncertainty as indicated in Figure ^ f° r 
the two high-mass bins. 

4. THE AGE DEPENDENCE OF THE CONDITIONAL 
LUMINOSITY FUNCTION 

In this section, we study the dependence of the occu- 
pation of galaxies on the halo formation time from the 
SAM and HYD galaxy formation models. The results are 
presented in terms of the CLF $(L|M), which is the ex- 
pected n umber of galaxies p er unit luminosity in a halo of 
mass M l|Yang et al.ll2003f) . Only galaxies brighter than 
Mr = —17.0 are considered in our study, where Mr is 
the i?-band absolute magnitude. 

The i?-band CLFs from the SAM and HYD models 
are shown in the left two columns of Figure [21 The drop 
below Mr ~ — 18 indicates that the CLFs are only com- 
plete for Mr < —18. For clarity, halos in each mass bin 
are only divided into three age bins with an equal number 
of halos in each age bin. A clear dependence of the CLF 
on halo age is found in both models, and the dependence 
becomes stronger at the faint end of the CLF. In gen- 
eral, there are more faint galaxies and fewer bright ones 
in younger halos. At the faint end, the youngest one- 
third of halos could host 10 times more galaxies than the 
oldest one-third. At the bright end, the trend of more 
bright galaxies in older halos is clear in the HYD model 



in the two high-mass bins, and it tends to disappear at 
low halo mass (~ 10 11 h~ 1 M G ). For the SAM model, the 
trend at the bright end is not apparent. The subtle dif- 
ference between the SAM and HYD models at the bright 
end could be due to the facts that the dust extinction 
effect is not included in the HYD model and that the gas 
cooling is switched off for halos of circular velocity larger 
than 350 kms -1 in the SAM model. 

In the central two columns of Figure |3 the CLFs are 
separated into contri butions from centr al and satellite 
galaxies, as is done in lZheng et alJ ((2005). At the bright 
end, the CLF is dominated by the bump caused by cen- 
tral galaxies. In older halos, central galaxies appear to 
be brighter, and the central galaxy bump as a whole 
shifts to higher luminosity. There are two reasons for 
central galaxies being brighter in older halos. First, in 
older halos, star formation is triggered at an earlier time 
and may last longer. Second, the accretion of satellite 
galaxies into older halos happens earlier, and there is 
enough time for them to merge onto central galaxies. 
Both scenarios increase the amount of stellar components 
in the central galaxies of older halos. To show that this 
is indeed the case, we plot the conditional stellar mass 
function (CSMF) *(M star |M) in the right two columns 
of Figure |21 separating into contributions from central 
and satellite galaxies. The CSMF is not complete at 
the low stellar mass end (M sta r < 10 9 /i -1 M Q ), which 
explains the steep drop at this end. It can be clearly 
seen that central galaxies in older halos have more stars. 
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Young stellar populations are more luminous than old 
ones, and this reduces, but not completely erases, the 
luminosity difference in central galaxies of young and old 
halos. As a result, the difference in the central galaxy 
CLFs in halos of different ages is not as large as that 
in the central galaxy CSMFs. In the lowest mass bin, 
the HYD model may have more stars formed recently, 
which broadens the central galaxy CLF in young halos 
and minimizes the difference among the central CLFs in 
halos of different ages. 

The faint end of the CLF, where the dependence on 
halo age becomes much stronger, is dominated by satel- 
lite galaxies. The correlation of the number of satel- 
lite galaxies with the halo formation time reflects the 
balance between the accretion and de struction processes 
ijGao et alJl2004HZentner et alJ feOOS). We find that, at 
a fixed mass, halos that formed earlier start to accrete 
satellites (subhalos) earlier. Then these satellites have 
a longer time to experience the orbital decay caused by 
the dynamical friction and will more likely merge with 
the main central halos. Subhalos surviving at z ~ are 
dominated by those accreted recently. Therefore, even 
though young halos accrete their mass later than old ha- 
los, they have much more surviving satellites. 

Since the stellar masses in central and satellite galaxies 
have opposite dependences on the halo formation time, 
we may expect that the total stellar mass in a halo of 
fixed mass is a weak function of the halo age. Defining 
the total stellar mass in a halo as the sum for galaxies 
with M st ar > 10 9 /i _1 Mq, we find that, in general, there 
are somewhat more stars in halos formed earlier, but 
the difference between the average stellar masses in the 
youngest and oldest halos is within a factor of 2. 

5. DISCUSSION AND CONCLUSIONS 

In this Letter, we investigate the dependence of the 
galaxy contents on the halo formation time using two 
galaxy formation models, one with a semi analytic 
method utilizing the realistic halo assembly history from 
A^-body simulations and the other with a smoothed par- 
ticle hydrodynamics simulation including radiative cool- 
ing, star formation, and energy feedback from galactic 
winds. Even though our techniques for building merger 
trees as well as our simulation methods are independent, 
we confirm the dependenc e of halo clusterin g on the for- 
mation time reported bv iGao et al.l l)2005[) . We study 
the CLF and CSMF of galaxies as a function of the halo 
formation time and find that galaxies inside halos have a 
strong correlation with the halo formation time. In halos 
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that form earlier, there are fewer satellite galaxies, and 
their central galaxies appear to be more luminous, which 
can be understood through the early onset of star for- 
mation and the merging of satellites onto central galax- 
ies. The dependence of the halo occupation number of 
satellite galaxies on the formation time in this study is 
simila r to that of subhalos presente d in iWechsler et all 
l)2005|) (also see lZentner et al.ll2005|) . which is not sur- 
prising since satellite galaxies reside in subhalos. 

The HOD or CLF model assumes that the statistical 
distribution of galaxies inside halos depends only on halo 
mass. However, we note that this assumption does not 
mean that at fixed halo mass, the occupation number is 
independent of halo age. The correlation between the 
occupation number and the halo formation time exists 
even if the halo asse mbly history is base d on the excur- 
sion set theory (e.g-. lZentner et aDl2005|) . and the distri- 
bution of the halo age at fixed mass could provide the 
scatter that goes into the probability distribution of the 
occupation number at fixed halo mass [e.g., P(N\M) in 
the HOD model] . It is the dependence of halo clustering 
on halo age that makes the above correlation potentially 
important in interpreting galaxy clustering data. Given 
the age-dependent halo clustering, the age dependence 
of HOD/CLF constitutes a key ingredient for a more ac- 
curate halo-based model of galaxy clustering. In view of 
this, a detailed study of the correlation between galaxy 
contents and halo formation time, like that performed in 
this Letter, is desired. 

Depending on the application, the current version of 
the HOD/CLF model without taking into account the 
non-Markovian nature of halo formation, can still give 
accurate descriptions fe.g.. lYoo et alJl2b05D . We reserve 
a detailed investigation for future work, to see to what 
extent it is adequate in different applications. 
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